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ABSTRACT 

The  theory  of  the  thin  shell  hollow  sphere  piezoceramic 
underwater  sound  radiator  is  investigated  in  detail.  Start¬ 
ing  with  the  equations  of  state  an  equation  of  motion  is 
derived  for  the  case  of  pure  radial  motion.  Formulas  are 
deducedfor  coefficient  of  electromechanical  coupling,  veloc¬ 
ity,  power  conversion,  mechanical  Q,  transmitting  pres¬ 
sure,  impedance,  receiving  response,  and  efficiency.  Design 
charts  based  on  the  properties  of  two  popular  ceramics 
(a  BaTi03  -  CaTi03  mix  and  a  PbTi03  -  PbZr03  mix)  are 
also  presented.  The  limits  of  validity  of  the  derived  equa¬ 
tion  are  also  discussed. 
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THEORY  OF  A  VERY  THIN  PIEZOCERAMIC  HOLLOW  SPHERE 
UNDERWATER  SOUND  RADIATOR 


INTRODUCTION 

The  vibrations  of  an  isotropic  elastic  sphere,  or  shell  bounded  by  two  spherical  sur¬ 
faces,  have  been  worked  out  in  detail  by  various  writers.  Poisson  (1)  in  his  memoir  of 
1839  discussed  the  radial  vibration  of  a  solid  sphere.  Jaerisch  (2),  Chree  (3),  and  Lamb 
(4)  investigated  problems  in  the  free  and  forced  vibrations  of  solid  and  hollow  spherical 
shells.  Bromwich  (5),  Jeans  (8),  and  Love  (7)  extended  these  investigations  to  include 
geophysical  phenomena.  Laura  (8)  worked  out  expressions  for  the  effect  of  a  liquid 
medium  on  a  vibrating  sphere.  In  more  recent  times  the  problem  of  the  acoustic  vibra¬ 
tions  of  a  polarized  electrostrictive  hollow  sphere  in  a  liquid  medium  under  electrical 
excitation  has  been  studied  by  Rosenthal  and  Baerwald  (9). 

The  following  analysis  recapitulates  the  work  of  Rosenthal  and  Baerwald  and  proceeds 
in  greater  detail  to  derive  results  not  found  in  their  report.  The  usefulness  of  many  of 
the  derived  equations  is  improved  by  presentation  of  design  graphs  in  Appendix  A. 


THE  STRUCTURE 

The  structure  (Fig.  1)  consists  of  a  hollow  sphere  of  polarized  electrostrictive 
ceramic  made  of  two  hemispheres  cemented  at  the  equator.  Each  surface  of  the  sphere 
is  completely  electroded  with  fired  silver  paste.  Polarization  of  the  ceramic  is  in  the 
direction  along  a  radial  line.  The  structure  is  driven  into  forced  vibration  by  the  appli¬ 
cation  of  an  alternating  electric  field  (E3)  applied  across  the  thickness  of  the  sphere. 


LIQUID  MEDIUM 
</>„  c„> 


Figure  1  -  Piezoceramic  hollow  sphere 
underwater  sound  radiator 
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The  vibrating  sphere  radiates  sound  into  an  infinite  medium  of  density  pm  and  sound 
velocity  (compressional  waves)  c„. 


EQUATIONS  OF  STATE 


The  analysis  will  be  conducted  in  the  orthogonal  spherical  system  of  coordinates,  e 
(colatitude  angle),  <p  (azimuth  angle),  and  r  (radius).  For  simplification  in  writing,  when 
these  appear  as  subscripts  they  will  be  designated  by  the  numerals  1,  2,  3  respectively.  It 
is  assumed  from  the  onset  that  the  thickness-to-exterior  radius  ratio  t/R„  is  considerably 
less  than  unity,  and  that  the  variation  of  mechanical  fields  (stress)  and  electrical  fields 
(electric  field  gradient)  with  the  radial  coordinate  is  negligible.  The  stress  field  and 
electric  field  are  thus  two  dimensional,  the  planar  character  of  which  will  be  emphasized 
by  placing  a  bar  over  all  symbols  representing  the  material  constants  (stiffness  constant 
5U ,  piezo  modulus  ?31,  etc.).* 


The  appropriate  set  of  equations  of  state  should  begin  with  the  physical  fact  that 
electromechanical  conversion  of  energy  is  obtained  by  exciting  a  polarized  electrostric- 
tive  material.  However,  it  is  more  convenient  to  treat  the  active  ceramic  as  an  artificial 
piezoelectric  crystal  and  assign  to  it  certain  "effective  moduli"  whose  form  and  magni¬ 
tude  are  determinable  by  conventional  piezoelectric  techniques.  In  accordance  with  this 
procedure  we  describe  the  planar  electromechanical  state  of  the  polarized  ceramic  by  a 
set  of  linear  piezoelectric  equations  in  the  form  suggested  by  Mason  (10).  These  equa¬ 
tions  state  that  the  planar  stresses  (Tj,  t2)  depend  upon  the  state  of  strain  (S„  S2)  and 
upon  the  electric  field  (e3)  in  an  infinitesimal  volume  of  material  anywhere  inside  the 
ceramic  in  the  following  way: 


Ti  =  5u  si  +  *12  S2  _  *31  E3 
tj  =  *12  si  +  *11  S2  "  531  E3 


(la) 


The  dependence  of  the  electric  displacement  vector  D3  on  the  electric  field  and  the  state 
of  strain  is 


03  -  c33  E3  +  #3i(Si  +  S2)  (lb) 

where  e3s3  is  the  dielectric  constant  at  zero  strain.  A  positive  strain  is  defined  as  one 
that  accompanies  a  state  of  mechanical  tension  in  a  body  free  to  move.  A  positive  elec¬ 
tric  field  vector  is  one  that  points  in  the  same  direction  as  the  polarizing  electric  field 
vector. 


The  mechanical  equations  (la)  and  the  electrical  equation  (lb)  are  Independent  state¬ 
ments  of  laws  of  nature  having  in  each  case  the  same  two  independent  variables  E  and  s. 
If  one  of  these  variables  (say  S)  is  eliminated  by  substitution,  a  single  equation  of  state 
is  obtained  which  describes  conditions  at  any  material  point  of  the  ceramic  due  to  both 
electrical  and  mechanical  fields.  Equation  (lb)  becomes 


°3  =  *33 


1  ♦ 


25 


3 


31 


*  *u>«sa 


E3  * 


*3l(T,  ♦  Ti) 


■11 


+  z 


12 


(2a) 


In  a  similar  way  the  variable  e3  may  be  eliminated  between  the  equations,  and  Eqs. 
(la)  become 


•A  list  of  symbols  appears  at  the  end  of  the  report. 
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T,  +  T,  =  (c 


_  E 
11 


_  E 


Clj) 


1  + 


2 

2*31 


(CE 

VC11 


5  E  w 
'12' 


33 


(S,  +  S4)  - 


2e 


31 


(2b) 


33 


It  is  noted  in  each  of  Eqs.  (2a)  and  (2b)  that  the  magnitude  of  the  material  constants 
e3s3,  a* ,  and  cfj  is  modified  by  the  same  group  of  assembled  symbols.  This  group  is 
definea  in  the  literature  (11)  as  the  material  static  mixed  planar  coefficient  of  electro¬ 
mechanical  coupling,  (k  2)  ;  that  is, 


2  2S32t 

<k,)- '  «,*,  *  *?,*& 


(3a) 


A  related  parameter  of  electromechanical  coupling  that  will  prove  of  use  in  later  discus' 
sions  is  (k  J)  ,  where 

Hon 


<kA.. 


<kp2)  , 

K  ais 


1  ♦ 


O  . 

K  Bis 


(3b) 


For  convenience  in  writing  we  can  reduce  the  number  of  terms  appearing  in  Eq.  (2b) 
by  recalling  that  ?1E1/5n  '  v  (Poisson's  ratio)  and  oy  defining  a  new  constant  CjD,  such  that 

5n  *  *iV>  +  fkn>.  J  • 

The  tangential  strain  S,  (St  =  s2)  for  a  purely  radial  symmetrical  mode  of  motion,  is 
simply  the  radial  displacement  j  divided  by  the  external  radius  Re .  Equation  (2b)  may 
therefore  be  written 


h  =  2c®  (1 


♦  v) 


S L  - 
R„ 


2e 


31 


33 


(4) 


In  this  equation  the  coupling  of  the  electrical  field  to  the  mechanical  is  implicitly 
contained  in  the  open  circuit  stiffness  constant  ctDt.  While  it  is  always  possible  to  use 
Eqs.  (la)  in  the  form  shown  as  an  aid  in  determining  the  radial  displacement  (,  the 
displacement  so  determined  would  be  a  first  approximation  only,  since  Eq.  (lb)  would 
still  remain  to  be  satisfied.  The  method  of  initial  substitution  used  here  is  convenient, 
though  superficial,  and  has  the  disadvantage  of  changing  the  field  variable  E3  to  the  elec¬ 
tric  displacement  variable  D3 .  However,  the  dynamic  equation  of  motion  is  directly 
deducible  from  Eq.  (4)  without  the  preliminary  inversion  of  factors  that  accompanies  the 
use  of  equations  involving  the  piezo  modulus  From  this  point  of  view  the  method 
leading  to  Eq.  (4)  is  preferable  to  other  possible  methods,  though  the  end  results  of  all 
correct  methods  will  differ  one  from  the  other  in  no  material  way. 


EQUATION  OF  MOTION  AND  SOLUTIONS 

The  stress  equation  of  dynamic  equilibrium  for  the  case  of  pure  radial  motion  in 
spherical  coordinates  is  (12) 


3T33  1 

3r  +  r 


3T 


31 


30 


[2T 

r  tin  0  30  r 


.J_ 


33  *  TU 


r22  +  *31  COt  01  ♦  pjr  =  Pt 


(5) 


All  the  terms  appearing  in  this  relation  are  forces  acting  upon  a  unit  of  (interior)  volume 
of  material  with  density  pc.  In  p  trtlcular,  the  quantity  Ff  is  the  body  force  per  unit  of 
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mass  arising  irom  some  external  agency.  This  differs  in  its  nature  from  the  elastic 
forces  accompanying  deformation,  which  are  described  by  the  stress  derivatives  above. 


Under  the  assumptions  of  thin  shell  planar  analysis  (restricted  to  pure  radial  motion) 
the  forces  stemming  from  the  stresses  t33,  t3 3,  and  T32  and  their  derivatives  are  to  be 
neglected.  In  addition  the  external  agency  of  (acoustic)  fluid  pressure  p  on  the  shell  of 
thickness  t  is  considered  to  induce  a  body  force  per  unit  of  volume  equal  in  magnitude  to 
P/t .  The  direction  of  this  "body  force"  must  be  carefully  chosen  to  agree  with  the  physi¬ 
cal  situation  of  a  pulsating  shell.  When  the  planar  stresses  Tn  and  t22  are  positive 
(tension),  the  radial  displacement  4  is  positive,  since  the  shell  is  expanding.  Under  the 
action  of  elastic  restoring  forces  the  acceleration  d2.j/dt2  must  however  be  negative, 
that  is,  toward  the  origin  of  coordinates.  Simultaneously  the  external  acoustic  pressure 
is  negative  (compression),  and  induces  a  negative  acceleration.  The  equation  of  motion 
therefore  for  a  very  thin  spherical  shell  of  external  radius  ro  reduces  to 


-rrn  ♦  t22) 


10) 


The  time  variation  of  4  has  not  been  specified,  and  while  pulsed  excitation  of  the 
sphere  would  prove  of  interest,  it  is  the  harmonic  response  that  is  of  principal  utility. 

In  addition,  since  we  have  made  acoustic  pressure,  p,  positive  when  it  applies  to  deficit 
pressure  (tension),  contrary  to  the  conventional  formulas  of  acoustic  theory  which  make 
excess  pressure  positive,  we  modify  the  latter  accordingly  and  write  the  known  specific 
acoustic  Impedance  of  a  pulsating  sphere  with  an  additional  negative  sign;  that  is,  we 
write  (13) 


-P.c.kr  (kr  +  j) 

P~kV 


v  , 


(7) 


where  k  is  the  wave  number,  and  v  is  the  medium  velocity  here  identified  with  j u£. 
Substituting  Eqs.  (4)  and  (7)  into  (6)  and  collecting  terms,  we  obtain 
,  *.cwkR0(kR0  +  j)  2*31 

^  +  td  .**o+  *n(1  ^  j*  -  TJ77  • 


(8) 


Equation  (8)  is  a  differential  relation  expressing  the  forced  excitation  of  the  sphere  by  an 
electromechanical  fore**  per  unit  volume  fo  equal  to  the  quantity  2*31  V'jjRo- 
radial  mechanical  displacement  Is  therefore 


£  1 

_ _ -.?■  — *  ...  -  - - 

/  u>  \  2  Pm  *0  /uA2 

1  '  \“n)  "  Pc  «  \%/ 

1  ~ 

■  '{$)* 

f 

V. 

L  '  "  '-j 

(») 


where 


(10) 


The  quantity  un  defined  by  Eq.  (10)  has  a  special  significance.  In  the  absence  of  an 
acoustic  load  (l.e.,  when  pm  -  0)  wn  is  that  frequency  which  produces  an  infinite  displace 
ment,  internal  dissipation  neglected.  It  is  therefore  the  ideal  mechanical  resonant  fre¬ 
quency  of  a  radially  pulsating  thin  shell  sphere  under  open  circuit  conditions. 
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A  noteworthy  feature  of  Eq.  (9)  is  that  the  displacement  is  a  complex  number,  result¬ 
ing  from  the  existence  of  velocity-dependent  and  acceleration-dependent  loads.  The 
absolute  magnitude  of  this  complex  displacement  is  found  to  be 

Is  I  --  ;s$'y* 


Further  deductions  from  Eq.  (8)  may  be  obtained  by  writing  it  in  a  more  revealing 
way.  Again  restricting  the  time  dependence  to  steady  state,  and  using  the  sinusoidal 
forcing  excitation  F„,  we  find  that 


r*j^n+  k*£0  =  Fa 

where 

■' =  '•  [“§  «.  '»’«.’>] 

r.  =  k2  *<? 

~  tn  +  kJRoJ> 

2c,D,n  +  v) 
k  =  — —  • 

O  2 
Ko 


(12a) 

(12b) 

(12c) 

(12d) 


Equation  (12a)  is  recognized  as  the  steady  state  differential  equation  of  forced  motion 
of  a  velocity-damped  mass-spring  system,  with  frequency-dependent  constants.  The 
solution  is  easily  seen  to  be 


F0  CO»  (hit  -  $) 

J r*  -  k* 


tan  <?■ 


r  oi 

-J*m*  *  k* 


(13a) 

(13b) 


The  mechanical  resonant  frequency  In  the  liquid  medium  is  defined  as  that  fre¬ 
quency  at  which  the  phase  angle  «  between  F„  and  (  is  */Z  That  Is, 


or 


7  R 


r — X 

i  ♦  —  n 

] 

{ 

1  ) 

J 

(14) 


The  difference  between  the  mechanical  resonant  frequency  in  air,  Eq.  (10),  and  the 
mechanical  resonant  frequency  in  a  liquid,  Eq.  (14),  thus  depends  upon  the  parameter  in 
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brackets.  This  latter  quantity,  which  may  be  thought  of  as  a  factor  by  which  the  density 
of  the  ceramic  is  to  be  increased,  is  itself  frequency  dependent,  and  hence  not  deter¬ 
minate,  unless  the  magnitude  of  »K  is  already  known.  Equation  (14)  is  therefore  best 
solved  by  trial  and  error.  However,  in  any  given  situation  it  will  be  convenient  to  assume 
o)R  *  a)n  in  order  to  compute  this  density  modifying  factor.  Repeated  trials  will  quickly 
yield  the  true  value  of  cjr. 

Further  deductions  from  Eq.  (13a)  may  be  made  by  proceeding  conventionally  to 
define  the  quality  factor  Q  *  and  the  resonance  factor  n  in  such  a  way  that  we  may  write 


F.  cos  (ut  -  <a) 


k‘  l"1 ' *4 


(17b) 


By  suitable  maximising  procedure*  the  frequencies  of  maximum  displacement  (w,) 
and  maximum  acceleration  )  are  derived  to  be 
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Because  of  the  simplicity  of  the  equation  of  motion,  Eq.  (12),  the  derived  expressions, 
Eqs.  (15),  (17),  and  (18),  were  quickly  obtained,  since  they  are  standard  results.  An  addi¬ 
tional  standard  result,  namely,  the  temporal  damping  factor  ad,  may  be  obtained  with 
similar  ease,  since 


and  thus 


ad 


r 

2m* 


P  c 


PA 


ko  R  ‘ 


!  *  kR2  R02  * 


H o 


(19) 


The  reduction  of  the  problem  of  the  radially  pulsating  sphere  to  that  of  a  system  with 
a  single  degree  of  freedom  is  a  convenient,  though  not  wholly  satisfactory,  approxima¬ 
tion.  Since  the  stiffness  per  unit  volume  is  k*,  the  total  stiffness  is  evidently  4nRoJ  tk*, 
a  quantity  proportional  to  the  thickness  t.  Having  idealized  the  elastic  system  to  the 
extent  of  allowing  t3/R0J  to  be  very  small,  it  is  seen  that  in  the  limit  the  hollow  sphere 
will  have  no  stiffness  at  all  (t  %  0).  It  will  resemble  a  membrane,  which,  for  any  given 
suifa^e  force,  yields  extremely  large  displacements.  The  acoustic  performance  of  such 
an  idealized  structure  will  be  extraordinary.  Due  caution  is  therefore  to  be  observed  in 
the  intended  uses  of  the  following  performance  formulas. 


ELECTRICAL  IMPEDANCE 


The  solution  of  the  equation  of  motion  has  led  to  an  expression  for  mechanical  dis¬ 
placement  i  in  terms  of  electric  displacement  D3  (Eq.  9).  Substituting  this  into  Eq.  (lb) 
and  solving  for  E3  we  obtain 


For  the  case  of  a  very  thin  shell,  the  current  I3  is  j ^4-n-Rj*  d3,  and  the  dependent 
voltage  Vj  is  approximately  E3t .  Hence  we  may  rewrite  Eq.  (20)  as  a  typical  impedance 
relation: 


h 


(21) 


where 


C*  *  « 


•t 


A  numerical  computation  using  Eq.  (21)  will  prove  instructive.  We  assume  a  hollow 
'  " - - - - - * - - - -  "  "  J'  *0,2S,*fn  c. 


shell  (In  water)  of  lead  titanate- lead  zirconate  ceramic  for  which  (k*) 
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at  air  resonance  is  4.0,  R0/t  -  10,  and  pw'pc  =  1/7.6.  If  the  drive  frequency  is  the  same 
as  the  mechanical  resonant  frequency  in  air  (if  u>  -  un)  then 

I3 

--•*  0.596  ja>C  +  0. 3S1  uiC  -  j  Bj  ♦  G,  . 

V3 

As  seen  from  the  electrical  terminals  the  sphere  appears  as  a  capacity  in  parallel 
with  a  resistance.  The  electrical  Q(  =  B,/G|Vis  approximately  1.5.  Due  to  the  inertial 
effect  of  the  water  the  susceptance  j is  reduced  to  0.6  of  its  low  frequency  value,  a 
remarkable  reduction,  stemming  from  the  high  coupling  factor. 


ACOUSTIC  POWER  DELIVERED 

From  Eq.  (12)  it  is  apparent  that  the  force  per  unit  volume  dissipated  in  the  resist¬ 
ance  r*  is  jr*.  Since  every  particle  of  the  sphere  moves  with  the  same  displacement, 
the  total  force  in  the  direction  of  the  velocity  4  is  4^Ro2  t  4  r*.  Recalling  that  D,  is  an 
independent  variable  (constant  current  source)  we  employ  Eq.  (17a)  and  write  an  expres¬ 
sion  for  the  maximum  acoustic  power  delivered  to  the  medium: 

=  (4  t R  2 1  )r*  4  £ 


16  v  Dj2  -  C) 

S  2  • 2  .  (  •  k*\ 

<33  r  +  H 


Since  the  applied  current  I3  is  equal  to  4*  Ro2  jwD3,  and  the  blocked  capacity  C*  is 
4ttro2  e  */t,  Eq.  (22)  can  be  rewritten  in  terms  of  I3: 

16rr  /  I3 \2  r*  sin2^t  -  Q) 


■’  *  («»•  -  $ 


It  is  frequently  desirable  to  have  an  expression  for  the  acoustic  power  in  terms  of 
the  applied  (or  source)  voltage  V3.  Regarding  the  source  current  I,  as  a  constant  source 
in  parallel  with  its  admittance  j^C*,  both  supplying  the  radiation  resistance  r*  with 
power,  we  see  that  the  equivalent  voltage  source  is  V3  in  series  with  the  impedance 
l/jctiC’,  both  supplying  the  same  resistance  r’ .  Hence  |l3/^C5 1  is  equal  to  V3,  the  peak 
applied  voltage. 

Employing  Eq.  (3b)  we  can  write  Eq.  (23)  in  terms  of  the  coupling  factor  (kp2)^  • 
877(k2)  cE(l+v)(l  +  (R2)  )  €,S,  r*V32  sin2(a't  -  £) 

T>  P  hom  *  * _ P  Mix  _ 

.  I  .  u-\’l  '  (2J> 

At.  the  frequency  of  velocity  resonance  (0  =  0)  the  rms  acoustic  power  into  the  water 
is  therefore 

4*(k2)  SE(1  *v)(l  +  (k02>  ]  e»(l  +  kR2  R  2)  V  2 

(P.)  = — - rjrr— - -•  (25) 

“mu*  V 


It  is  to  be  remembered  that  V3  in  Eqs.  (24)  and  (25)  is  the  peak  voltage. 
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Equation  (25)  shows  that  the  acoustic  power  output  at  velocity  resonance  is  approxi¬ 
mately  independent  of  the  radius  R0  To  reveal  this  fact  we  note  that 


L.J 


R2  =/«y 
0  \c«/  \%j 


2  (1  ♦  u) 


When  -j  =  <x, ,  we  find  that 


kR  RoJ  =  2  CjEj  (1**)  Vc  cw2 


which  is  a  factor  independent  of  R0.  When  R  <<  n  this  conclusion  fails. 

While  the  total  resonant  power  per  voltage  squared  is  approximately  independent  of 
ROI  the  surface  power  density  (watts/cm2)  is  not.  Hence  a  large  resonant  sphere  (a  low- 
frequency  sphere)  cat)  radiate  much  more  energy  than  a  small  resonant  sphere,  because 
it  can  tolerate  a  larger  voltage  (that  is,  for  equal  maximum  surface  power  density,  its 
voltage  drive  will  be  very  much  larger). 

An  illustrative  computation  using  Eq.  (25)  will  prove  of  interest.  We  select  a  lead 

titanate-lead  iirconate  material  for  which  fk  2 ) h  -  0.23,  cjEj  r  91.3  x  10\  v  =  0.3, 

k  '*)  ,  t s  -  «■ 

P  ha.  33 

2  '  01.3  x  109a.3) 


- - - - - - —  - >  -  --  -  » 

1  +  (k  2)  s;  l  +  (k  2)  ,  i  s  =  j.hq  .  io-o,  and  oc  -  7. fi  x  103.  At  velocity  resonance 


k  2  R  2 

kR  ko 


=  13.0  . 


Hence 


7.6  -  103  *  1.552  x  io6 

P  -  477  x  0.23  x  91.3  X  1Q9(1.3)  (1.23)  7.89  x  IQ-9  x  j4 
’  "  1.55  x  106  >■•  13 

-  2310  «  in-6  V2  (rms  power:  V  =  peak  volts)  . 

For  a  sphere  5  in.  in  diameter  and  1/4  in.  thick  driven  with  850  peak  volts,  the  rms 
power  would  be  1670  watts  at  100-percent  electromechanical  efficiency. 

We  may  also  deduce  a  formula  for  the  voltage  required  to  achieve  a  surface  power 
density  i>  (watts/m2)  based  upon  the  restricted  case  of  Eq.  (25).  At  velocity  resonance, 
we  find 


v3  =• 


STb  R 


oo 

y  ho 


.  /[>  * 


fkp2) 


J~  E33  Cw(*  +  V  R02) 


(26) 


Example:  To  achieve  a  density  of  3  watts/cm2  using  the  5-in.  sphere  of  the  material 
noted  above  we  require 


v3  = 


i/Tx  3  x  io*  x  6,35  X  10-2 


0.48  / (1.23)(7.6)  7.89  x  io*9  x  1,5  x  103(14) 


=  825  peak  volts. 


ACOUSTIC  PRESSURE  IN  THE  FAR  FIELD 

The  expression  for  acoustic  power  obtained  in  the  previous  section  may  be  written 
symbolically  as  the  product  of  the  particle  pressure  and  the  conjugate  particle  volume 
velocity.  If  the  particle  velocity  Is  u,  the  latter  quantity  is  4wrV,  since  all  parts  of  the 
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sphere  vibrate  in  phase  with  each  other.  Now  in  the  far  field  at  distance  R  the  particle 
pressure  and  velocity  are  substantially  in  phase  and  are  related  to  one  another  through 
the  characteristic  impedance  of  the  medium  pwcw.  The  conjugate  volume  velocity  is  then 
4 11  R2p */£, <v  The  magnitude  of  acoustic  power  in  the  far  field  therefore  becomes 

ip  i  I  p  i 2 . 

*  x?  c 

■  w  » 

Equating  this  magnitude  of  acoustic  power  with  the  magnitude  in  Eq.  (24),  and  solving 
for  the  peak  acoustic  pressure  in  the  far  field,  we  obtain 


IP  I  =- 


71  <kp>  V3 

^  hoa 

I  CK  C1  +  ^ 

1  +  (k2)  .kV* 

L.  K  ipi  aj 

R 

/ 

j 

f  t 

'*2  +l 

f  .  k*\2l 

(  "/  1 
J 

(2?) 


At  velocity  resonance,  -  k*  4>,  and,  after  substituting  Eq.  (12c),  the  far  field 
magnitude  of  acoustic  pressure  reduces  to 


IP  I  =  - 


^  0te)ho.V3  1^)11  ♦  «3Vl  +  K2  V> 


Wo 


(28a) 


ipi  - 


/T  (kp\0.r3  ^ (l  +V)il  +  ‘aV1  +  kR  ^ 


R 


kRRo 


(28b) 


At  velocity  resonance,  the  magnitude  of  the  far  field  pressure  is  approximately 
independent  of  the  radius  of  the  sphere,  provided  the  frequency  in  the  liquid  medium  does 
not  differ  appreciably  from  the  frequency  in  air.  This  will  be  the  case  for  radius-to- 
thickness  ratios  (ao/t)  less  than  20  A  sample  calculation  using  the  constants  of  the  5-in. 
lead  titanate-lead  zirconate  sphere  mentioned  above  (in  which  R0/t  =  10)  shows  the  value 
of  the  quantity  |p|/V3  to  be 

[p |  _  Jj  x  0.48  /91.3  x  iq»  x  1,3  x  1.23  «  7. B9  *  IQ-9  *  14 
V3  ~  R  yT3 


23. 8  newt  on  _  47.6  do  rei.  1  uoar  at  1  meter 
R  meter-volt  ~  volt 


AI30,  since 


<jC*  =  2v  *  14  x  103  x  57,916  *  10' 12 


=  465  *  10*5 


we  have 

!p|  5100  newton  94.2  db  ref.  1  Mbar  at  1  meter 

I  R  meter  ampere  ampere 


These  are  the  maximum  values  to  be  expected  since  the  efficiency  of  electromechanical 
conversion  is  assumed  to  be  100  percent. 
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RECEIVING  RESPONSE 

As  long  as  the  square  of  the  acoustic  pressure  in  the  far  field  is  strictly  propor¬ 
tional  to  applied  electrical  power  {linear  operation)  the  acoustic  performance  of  the 
sphere  will  obey  the  reciprocity  law  of  linear  transducers.  Defining  the  reciprocity 
factor  J  for  distance  ft  and  wavelength  ft  as 


J  »  (]'  *  ~~  for  R  =  1  meter') 

J  ft.c.  \J  Pwc*  ) 


and  recalling  that  the  open  circuit  receiving  response  M  is  related  to  the  short  circuit 
transmitting  response  S  through  J,  we  can  write 

-  M  -  TS  -  ----  ^ 

fol  “  ”  J  ^  ~  c  I 

r  w  *3 

Replacing  V*  by  wt  r  k*  because  of  electromechanical  coupling,  we  intro¬ 

duce  Eq.  27  into  this  equation  and  obtain  the  open  circuit  receiving  response  of  the  sphere: 

|pt  ~  uc *  r  /  .t\Jl  ’  (29) 

I  *|/a+r  -U 

Numerical  Example:  At  14  kc  the  factor  J  has  a  value 

2ft  _  Jrr_ _ 4rr _ 

J  "  P„c„  _  a pm  =  2w  «  14  *  103  x  io3 
=  j  x  JO'6  . 

From  the  results  of  the  previous  computations  on  a  5 -in  sphere 

=  i  x  10*6  x  51J0  =  7.3  x  io-4  ,  7.3  x  io-3  *2*1 

Ip  I  7  newton  Mbar 

=  -82.7  db  ref.  1  volt/^bar. 


TRANSMITTING  EFFICIENCY 

The  dissipation  factor  r*  is  the  useful  radiation  load  delivered  through  the  agency  of 
the  applied  electrical  power.  If  we  now  assume  an  internal  mechanical  (parasitic) 
resistance  rm  proportional  to  velocity,  we  can  write  the  total  mechanical  dissipation 
factor  as  r*(l  +  x),  where  v  =  rm/r*.  The  rms  power  Plfm  absorbed  by  both  radiation 
and  mechanical  damping  is  found  by  replacing  r*  by  r*(l  +  v)  in  Eq.  (24),  and  taking  1/2 
the  resultant  magnitude  since  peak  voltages  are  employed.  The  rms  power  Pd  absorbed 
by  the  dielectric  resistance  of  internal  di  mping  b  is  given  by 

y  2  y  2 

Pd  =2T  t8n  (30) 


Adding  the  power  absorbed  mechanically  to  the  power  absorbed  electrically  we  obtain  the 
total  absorbed  power  PT.  The  efficiency  v  is  therefore 
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Appendix  A 
DESIGN  GRAPHS 


A  series  of  design  graphs  based  upon  the  formulas  derived  in  the  main  body  of  the 
report  and  on  the  material  constants  listed  in  Table  A1  is  presented  in  the  following 
pages.  Due  to  the  limitations  of  planar  analysis  as  applied  to  a  sphere  and  the  uncer¬ 
tainty  of  the  values  of  material  constants  in  a  given  (physically  real)  sphere,  it  is  best  to 
regard  the  numerical  results  found  from  these  graphs  as  convenient  guides  only,  not 
assuming  an  accuracy  better  than  ±7  percent.  A  few  comments  on  the  construction  of 
these  graphs  will  be  appropriate  here.  In  the  course  of  these  comments  we  will  consider 
a  typical  design  problem  illustrating  their  use. 


TABLE  A1 

Material  Constants  of  Polarized  Ceramics 


Ceramic  B 
(BaTiCK  -  CaTi03) 

PZT-4 

(PbTiO,  -  PbZr03) 

PZT-5 

(PbTi03  -  PbZr03) 

CjFj  (newton/meter2) 

127.9  x  109 

91.3  x  109 

74.7  x  109 

cf2  (newton/meter2) 

38.3  x  109 

27.9  x  io9 

23.0  x  io9 

S®  (meter  2/newton) 

8.62  x  io- 12 

12.05  x  io* 12 

12.05  ><  10'12 

<vh0. 

0.33 

0.48 

0.54 

d31  (coulomb/newton) 

-58  x  IO’12 

-97  *  io'12 

1 

h-» 

o 

X 

o 

K> 

e31  (planar) 

(newton/meter-volt) 

-9.45 

-10.38 

-12.01 

s3s3  (coulomb/meter-volt) 

9.64  x  10-’ 

7.89  x  10'9 

10.24  x  IO'9 

«3t3  (coulomb/meter-volt) 

10,7  x  l0-q 

9.7  -  IO*9 

13.22  x  IO'9 

pe  (kilogram/meter3) 

5.4  x  103 

7.6  x  io3 

7.6  x  IQ3 

■ 

(meter/ second) 

8.28  *  103 

6.2  *  103 

5.73  *  103 

0.3 

0.3 

0.3 

Figures  A1  -  Resonant  Frequency 

The  curves  in  Figs.  A1  are  numerical  values  of  Eq.  (10),  except  that  is  used  in 
place  of  c®.  The  resonant  frequencies  shown  are  correct  as  long  as  the  shell  is  very 
thin,  that  is,  as  long  as  Re/t  is  greater  than  10.  Let  us  begin  our  typical  design  problem 
by  requiring  a  PZT-4  sphere  (see  Table  Al)  to  vibrate  resonantly  in  water  at  1  kc.  From 
Fig.  A1  we  find  that  ro  is  35  Inches  in  air.  The  final  design  radius  will  be  smaller  because 
of  water  loading. 
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Figure  A2  -  Mechanical  Q 

Figure  A2  holds  for  both  materials  discussed  in  the  report.  It  is  based  upon  Eq.  (16b), 
for  which  the  condition  k2Ro2  »  1  holds.  In  our  problem  we  specify  a  mechanical  Q  of  2. 
From  Fig.  A2  we  see  that  R0/t  must  be  about  15,  and  hence  that  the  thickness  t  must  be 
2-1/3  inches. 


Figures  A3  -  Radial  Displacement 

Figures  A3  are  based  upon  Eq.  (11).  The  velocity  of  sound  in  water  is  taken  at 
cw  =  1.55  *  io3  meter/sec.  For  a  value  of  Rc/t  =  15,  we  judge  the  radial  displacement 
to  be  2.75  «  10  •9  meter/volt  at  u  =  <^n,  and  somewhat  higher,  say  2.90  *  10'9  meter/sec 
at  true  water  resonance.  The  surface  velocity  per  volt  is  therefore 

«  displacement  =  1v  x  103  *  2.75  *  10"9  =  1.73  *  10’*  m/volt-aec. 


Figure  A4  -  Surface  Area 

Figure  A4  gives  the  surface  area  of  spheres  as  a  function  of  diameter.  The  surface 
area  of  a  70-inch-diameter  sphere  is  99,000  cm2. 


Figures  A5  -  Power  Density 

The  curves  in  Figs.  A5  are  numerical  plots  of  Eq.  (26).  If  we  allow  a  surface  power 
density  of  1  watt/cm2,  we  will  need  approximately  6900  volts  (peak)  electrical  drive. 
Since  the  thickness  is  2-1/3  inches,  or  6.45  cm,  the  maximum  electrical  gradient  will  be 
818  volts/cm.  The  material  of  the  sphere  is  evidently  far  from  its  power  absorption 
potential,  since  a  value  of  2400  peak  volts/cm  could  be  used  with  impunity.  The  total 
power  for  100-percent  efficiency  and  for  the  coupling  cited  will  be  about  100,000  watts. 


Figures  A6  -  Power  Curves 

Figures  A6  are  based  upon  Eq.  (24).  We  have  plotted,  however,  the  rms  values  (1/2 
of  the  formula)  instead  of  the  peak  values.  In  our  design  problem,  let  us  require  the 
power  (for  100-percent  efficiency)  at  1/2  the  design  frequency,  that  is,  at  u>  /un  =0.5.  We 
find,  from  Fig.  A6b,  that  the  real  power  output  is  approximately  250/2285  or  0.09  of  its 
value  at  resonance.  Hence,  the  power  at  500  cps  will  be  about  11  kw. 


Figures  A7  -  Resonant  Pressure  Amplitude 

Figures  A7  are  based  on  Eqs.  (28a)  and  (28b).  From  them  we  find  that  for  a  drive  of 
6900  peak  volts  (76.8  db  ref.  1  volt)  the  source  level  at  1  meter  will  be  125  db  ref.  1  >ib«r 
(Fig.  A7A)  and  the  transmitting  response  will  be  +89.2  db  ref.  1  »bar  per  ampere  (Fig. 
A7b). 


Figure  A8  -  Receiving  Response 

The  curves  in  Fig.  A8  are  based  on  Eq.  (29),  computed  at  velocity  resonance.  For  a 
thicknes  of  2-1/3  inches  the  resonant  receiving  response  will  be  -65.0  db  ref.  1  volt/*jb*r. 
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Figure  A9  -  Blocked  Reactance  at  Resonance  (xB) 

The  curves  in  Fig.  A9  are  computed  from  the  simplified  assumption  of  XB  =  1/og  c\ 

At  velocity  resonance  uKc*  is  a  function  of  t/R0  only.  For  t/R0  =  1/lS  =  0.067,  XB  =  l20ohms. 


Figure  A10  -  Electrical  Qg  Curves 

For  a  given  material  the  coupling  factor  k  is  known.  Using  the  simplified  relation 
that  Qg  =  [1  -  (k  *''ho.]/(kpJ)hoI,  the  curves  for  the  two  materials  in  question  were  com 
puted  as  shown  in  Fig.  A10.  In  the  illustrative  example  Qu  =  2;  hence  QE  =  1.7.  Now 
assuming  that  qe  is  roughly  X^/R,  we  compute  R  *  70  ohms.  The  electrical  impedance  of 
the  sphere  is  therefore  approximately  70  -  j  120  at  velocity  resonance. 
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Figure  Ala  -  Resonant  frequency  in  air  of  thin  ceramic  spheres  versus  the  radius  (R„)  of  the  spheres 
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Figure  Ai.  -  Approximate  mechanical  Q  of  water  loaded  piezoelectric 
ce.amic  spheres  versus  tne  radius-to-thickness  ratio  (1?  A  )  This 
curve  applies  to  all  three  materials.  0 
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Figure  A3a  -  Radial  displacement  per  volt  versus 
frequency  ratio  «A>n  and  the  radius -to -thickness 
ratio  of  piezoelectric  ceramic  spheres 
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Figure  A4  -  Surface  area  of  a  sphere  versus  diameter 
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RESONANT  ACOUSTIC  SURFACE  POWER  DENSITY  8  (WATTS/CM* ) 


Figure  A5a  -  Peak  voltage  required  (100%  efficient:’')  versus  the  rms  value  of  the 
desired  acoustic  surface  power  density.  The  values  shown  arefor  PZT-4  and  PZT-5; 
Ceramic  B  values  can  be  obtained  by  multiplying  Vj  by  1.165. 
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RESONANT  ACOUSTIC  SURFACE  POWER  DENSITY  8  (WATTS/CM*) 


Figure  A5b  -  Peak  voltage  required  (100%  efficiency)  versus  the  rmi  value  of  A,  the 
desired  acoustic  surface  power  density.  The  values  shown  are  for  PZT-4  and  PZT-5; 
Ceramic  B  values  can  be  obtained  by  multiplying  Vj  by  1.165. 
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Figure  A6a  -  Reel  acoustic  power  (rmi  watts)  per  squared  peak  applied 
voltage  veraua  frequency  ratio  for  pieaoceramic  apheres  with  various 
R  /t  values 
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Figure  A6b  -  Real  acoustic  power  {rmi  *»Hi|  per  squared  peak  applied  voltage 
versus  frequency  ratio  -•%  for  piescceramic  spheres  with  various  R#''t  values 
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RESONANT  PRESSURE  AMPLITUDE  (db  REF.  BAR  AT  I  METER) 
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Figure  AT*  -  Reaonant  preaiure  amplitude  ver»u«  peak  applied 
voltage  (!00%  efficiency) 
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Figure  A7b  -  Reaonant  preteure  amplitude 
venue  the  radiue -to-thickm  ■■  ratio 
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Figure  A8  -  Resonant  open  circuit  receiving  response 
versus  shell  thickness  t 


LIST  OF  SYMBOLS 


Planar  stiffness  coefficients  (newtnn/mstar2) 

Capacity  at  zero  strain  (farad) 

Cerapresaional  velocity  of  30und  in  water  (meter/second) 

Piezo  modulus 

Electric  displacement  (coulomb/meter2) 

Outside  diameter  of  sphere  (meter) 

Superscripts  indicating  at  constant  electric  displacement  and  constant 
field  respectively 

Electric  field  (volt/meter) 

Planar  piezo  modulus  (newton/meter-volt) 

Body  force  per  unit  of  mass 

2*31°3/€33Ro 

Applied  current  (ampere) 

Reciprocity  factor  (meter  s/newton- second) 

Wave  number  (meter*1) 

Defined  by  Eq.  (12d) 

Wave  number  at  velocity  resonance  in  air  (meter*1) 

Mixed  coefficient  of  electromechanical  coupling 
Homogeneous  coefficient  of  electromechanical  coupling 
Wave  number  at  velocity  resonance  in  water  (meter  **) 

Defined  by  Eq.  (12b) 

Alternating  acoustic  pressure  (newton/meter2) 

Acoustic  power  (watt) 

Quality  Factor  (  =  mechanical  at  velocity  resonance) 
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Radial  coordinate 
Defined  by  Eq,  (?,2c) 

Internal  mechanical  resistance 
Distance  in  far  field  (meter) 

Outside  radius  of  sphere  (meter) 

Superscript 

Planar  strains  (meter/meter) 

Compliance  (meter  2/newton) 

Shell  thickness  (meter) 

Planar  stresses  (newton/meter2) 

Stresses  (newton/ meter2) 

Radial  velocity  in  steady  state  (meter/second) 

Applied  voltage  (volt) 

Open  circuit  voltage  (volt) 

Temporal  damping  constant  (Eq.  19) 

Dielectric  constant  at  zero  strain  (coulomb/meter-volt) 

Dielectric  constant  at  zero  stress  (coulomb/meter-volt) 

Spherical  coordinate 

Poisson's  ratio 

Radial  displacement  (meter) 

Radial  velocity  (meter/second) 

Density  of  ceramic  (newton- second 2/Wter4) 

Density  of  water  (newton- second  Vmeter4) 

Spherical  coordinate  and  phase  angle 

r«A* 


oi  Frequency  (radian/second) 

wn  Resonant  frequency  in  air  (radian/second) 

oiR  Resonant  frequency  in  water  (radian/ second) 


